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Abstract 

A system of a metastable phase with several sorts of the heterogeneous 
centers is considered. An analytical theory for the process of condensa- 
tion in such a system is constructed in dynamic conditions. The free 
energy of formation of the critical embryo is assumed to be known in the 
macroscopic approach as well as the energy of solvatation. The recurrent 
procedure of the establishing of the characteristic times and lengths is 
presented. The semiempirical method of the description of the period of 
essential formation of new phase embryos on the heterogeneous centers of 
some sort is given. 
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1 Introduction 



The theory considered here wiU be based on the capillary approximation of the 
height of the activation barrier. This publication can be regarded as the direct 
continuation of the publications and is going to fulfil the program 

announced in . All necessary bibliographical cites can be found in jl^ . 

It is necessary to stress that in the nature the more spread conditions arc the 
external conditions of the dynamic type. In such conditions the external action 
on the volume of the system, the temperature and the pressure are changed 
in some continuous smooth way. The process of condensation violates these 
characteristics. But this process itself takes place under the action of external 
influence on the metastability of the system. These conditions were considered 
at first by Raiser j?). The theory of the heterogeneous condensation on the 
centers of one sort was done in |^ . 

Practically in the condensing system there are several types of the heteroge- 
neous centers of the different nature (see ]lO|). 

The theory when the processes of the homogeneous and heterogeneous nu- 
cleation are taken into account in the dynamic conditions was constructed in 
[|j . The problem of construction of the theory for condensation on the hetero- 
geneous centers of several sorts still remains actual. Here such a theory will be 
constructed. 

Here and in the further considerations all energy-like values are expressed 
in the units of the thermal energy kjjT (kf, is Bolzman constant, T is abso- 
lute temperature). The theory presented here will be valid also for the pure 
homogeneous process. 

We shall use the following physical assumptions analogous to jl^ , |jll| : 

• the thermodynamic description of the critical embryo, 

• the random homogeneous space distribution of the heterogeneous centers, 

• the free-molecular regime of the droplets growth, 

• the homogeneous external conditions for the temperature and for the pres- 
sure, 

• rather a high activation barrier of condensation, 

• the absence of the thermal effects. 

As far as the most interesting characteristics of this process are the numbers of 
the heterogeneously formed droplets of the different types we shall estimate the 
accuracy of the theory by the error of the obtained solutions for these values. 
The unit volume is considered. 

We assume the total number of the heterogeneous centers to be constant in 
time. 
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The indexes below the values may refer to some characteristic 

times and t**, correspondingly. 

The theory of the process appear to be rather long and complex. So we 
shall use some already known constructions |^ , Q as some blocks in this theory. 
Recall that the heterogeneous condensation under the action of the linear source 
of the metastability was considered in |^ . The same problem for the source of 
the "square" form was studied in ^. 

Suppose that in the system there are several sorts of the heterogeneous cen- 
ters distributed homogeneously in the space volume. We shall mark the total 
number of the heterogeneous centers by rjtot i where i corresponds to some sort 
of heterogeneous centers. The real values of the free heterogeneous centers which 
may be solvatated but aren't occupied by the super-critical embryos are marked 
by rji. The index i or j below the value marks the sort of the heterogeneous 
center. The absence of this index points that the formula is valid for an ar- 
bitrary sort of the heterogeneous centers. The density of the molecules in the 
equilibrium vapor is marked by rioo, the density of the molecules in the real 
vapor in the system is marked by n. The power of the metastability will be 
characterized by the value of the supersaturation 



We shall define the super-critical embryos as the " droplets" . Every droplet is 
described by the number of the molecules inside the droplet i' , or by the linear 
size p = v^l^ . Due to the free-molecular regime of the vapor consumption we 
have 

. -1 

where a is the condensation coefficient and r is some characteristic time between 
collisions of the given molecule of the vapor with the molecules obtained from 
the gas kinetic theory. 

Let us introduce some size z according to 



^ C^ar-^dt' (1) 



where is some characteristic time. The choice of is more complex than in 
[||. For every sort of the heterogeneous centers there will be some special i. 
It is chosen as the moment until which the half of the droplets formed on the 
centers of the given sort is already appeared. In such a case the variable z will 
be marked as Zi. 

Until the coalescence [||, which isn't considered here equation (|l|) ensures 
the growth of z in time and can be inverted 

tiz) = a-^f^^+U (2) 
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Hence, all values dependent on time become the values dependent on z and the 
relative size x = z — p can be introduced. For Zi the values introduced 
by the same procedure. 

During the whole evolution the droplet has one and the same value of the 
variable Xi. Considering t{x) (or ti{x)) as the moment when the droplet with 
the given x has been formed (as a droplet) we can consider all functions of time 
as the functions of x . Hence, we can see that the kinetic equation is reduced 
to the fact that every droplet keeps the constant value of x. To reconstruct the 
picture of the evolution one must establish the dependencies t{z) and C{x). 



2 The system of the equations of condensation 

To describe the action of the external conditions we shall introduce the value of 
the ideal supersaturation 

Qd = 1 

Due to solvatation of the heterogeneous centers the ideal supersaturation 
falls to the value 

$ = Cid - '""^ ' (3) 

Tloo 

where i is the number of the molecules of the condensated phase in the 
equilibrium heterogeneous embryo. 

One can assume that Ve i can be taken at ^ = C*- 

One can put approximately 

^ = Cid 

with a small relative error. Wc shall call $ the ideal supersaturation also. 
To construct the mathematical model we must formulate some statements: 

• (1) The main role in the vapor consumption (when this consumption is re- 
ally essential) during the evolution is played by the super-critical embryos, 
i.e. by the droplets. 

• (2) The quasistationary approximation for the nucleation rate can be ac- 
cepted during every period of the essential formation of the droplets on 
every sort of the heterogeneous centers (except the situations when the 
considered sort of the heterogeneous centers is exhausted). When this 
approximation is essential it is valid. Namely, it is not valid only in the 
situation when all the nucleus of the given sort become the centers of 
the droplets. We shall call this situation as the "exhausted centers sit- 
uation" (ECS). In this very situation the quasistationary approximation 
isn't necessary as far as the result of the process is obvious. 
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• (3) It is possible to change in the formaUsm of the iteration procedure 
in Q to the choice of t** which corresponds to the maximum of the su- 
persaturation. Really, the "correct" choice as the time when the half 
of the total number of the droplets is already formed begins to deviate 
essentially from the choice of t,, only in ECS when the result is already 
known. Certainly, the parameter in the linearization must be reconsidered 
if it is necessary]^. 

The justification of the second statement uses the estimate for the times tf 
of the establishing of the stationary state in the near-critical region which can 
be found in [|], Q (for the heterogeneous barrier the consideration is quite the 
same). Here it is necessary to remark that, certainly, there may exists some 
rather huge times tf. They correspond to rather big values of the number of 
the molecules in the critical embryos. The halfwidth of the near-critical region 
[|] can be estimated by the homogeneous value at the supersaturation corre- 
sponding to same value of the volume of the critical embryo. It is proportional 

2/3 

to the value near lyj^ . Here and in the further considerations the lower index 
"c" marks the values for the critical embryos. As far as the absorption ability 

2 /3 

is proportional to i/^ ■ and the size of the near-critical region is proportional 

2 /3 2 /3 

to v^'^ the value of tf is proportional to . Practically, the big value of 
the activation barrier AFi — Fi{vc) (in the case when Vc i is greater than the 
characteristic length Ax of the size spectrum) means that these kinds of the 
heterogeneous centers are excluded from the kinetic process. 

For the majority of the types of the heterogeneous centers the following 
approximations for the nucleation rate Ji are valid during the period of the 
essential formation of the droplets on the centers of the given sort 

J^ = MVtot ^, C*) exp(r, ~ )^ (4) 

C,* Vtot i 

where 

_ d\F{C) 
' ~ ~ — dC, — ^ ' 

The validity of these approximations is justified for the heterogeneous embryos 
with the monotonous interaction between the nucleus and the molecules of the 
condensated phase which is weaker or equal than the reciprocal to the space 
distance. 

Let /* i be the amplitude value of the size distribution of the heterogeneously 
formed droplets measured in the units of rioo- As far as the supersaturation 
and the number of the heterogeneous centers rjtot i are the already known values 
the stationary rate of nucleation J and the stationary distribution / can be easily 

^Here the simplest explanation is given. The precise analytical analysis which confirm this 
statement is also available. 
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calculated by some formulas derived in 



= / = „ -1/2 A ..A ■.^^T ?^ (6) 



where is the number of the molecules absorbed by the embryo in the unit 
of time, AeZ/ is the width of the equilibrium distribution 

^e^V= exp(-i^(z/)) 

and Ac iV is the halfwidth of the near-critical region 

21/2 



d2p ,1/2 
0^2 \u=W^ 



We shall mark by rioo^i the total number of the vapor molecules in the 
droplets formed on the centers of the sort " i" . To simplify the formulas we shall 
use 



Vtot i 

Let us choose for the simplicity the concrete conditions of the constant vol- 
ume and pressure of the system. This allows to cancel some unimportant pa- 
rameters and to exclude them from formulas. 

Using the conservation laws for the heterogeneous centers and for the molecules 
of the condensated substance we obtain for g^, 9i the following equations 



gi = f*i I {zi-Xif ex^{Ti^^^ — ^^)exp{-^i^r-^)Oidxi = Gi(^gj,0i) 

J — oo i C* i j 

(7) 

0, = exp(-/,,^ r exp(r,iM^^)exp(-r,^)dx,)^ 5,(5^5,) 

'fjtot i J — oo S* j S* i J 

(8) 

where /* i = Ji{ritot i,C* i)'^/C* iauoc,- 

These equations form the closed system of the equations of the condensation 
kinetics. This system will be the subject of our investigation. 

As far as we measure the accuracy of the theory in the terms of the error in 
the number of the droplets we define these values as the following ones: 

= j^tot ^{l - e,{z)) = Q,{e,) (9) 

The size spectrum can be found as the following one 

/, = , exp(-r,^) eM^^^^^^)8^ (10) 
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Note that system - (||) is rather more complex than in the case of decay. 
The specific feature is not only the presence of the factor exp(ri ' ) but 

also the necessity of the introduction of the separate scales for z and x. It means 
that the periods of the essential formation of the droplets on the centers of a 
different nature don't coincide or have any common regions. Such a situation 
can be guaranteed by the following inequalities 



f_ 

or 



< 1 (11) 
> 1 (12) 



/. 

for all arbitrary pairs of the sorts of the centers i and j. The situation 

^ ^ 1 (13) 

J*j 

will be at first excluded from consideration. Later the generalization will be 
presented. 

Note that inequalities ( |ll| ) - ( |l3| ) are rather unsensitive to the change of the 
supersaturation. Note that 

/, exp(-Fc + Fe) ^ exp(-Fc hom + Fe) 

where the free energy F is counted in the "absolute" scale with the zero point 
of the energy corresponding to the standard homogeneous zero point Fhom{v = 
1) = and Fj e are the energies of solvatation, So, we can write 

/* i ^ Vi tot exp(+Fie) , , 

f*, Vj tot exp{+Fje) 

Hence, we can see that when the free energy of solvatation is unsensitive to the 
supersaturation (as it occurs really in the nature) then the r.h.s. of the last 
equation doesn't depend on the supersaturation. 

So, the set of the heterogeneous centers {rji tot} initiates the set of the char- 
acteristic times and the characteristic supersaturations {Ci*}- Suppose 
that Q , > i if * > j- 



3 External supersaturation 

The analysis of the last section simplifies the treatment, but isn't absolutely 
necessary. The more rigorous recipe is to solve the separate problems 

$(x) ' 
{zi-Xi)^ exp{Ti —)exp{-Ti—^)0idxi = Gi{gi,9i) (15) 

-oo <,* i S>* i 



7 



Vtot i J ~oo S* i S* i 

This solution can be obtained by the methods from |^] and we needn't to 
describe it here. As a result we get the set {C* i}. The solution is valid only 
for I = 1 (C* 1 < C* AH other solutions are got only to construct the set 

{C* i\- Now they are useless. 

The further condensation can't be described by system ( |l5| ) - (16) because 
the droplets formed on the centers of the first sort consumes the vapor. This 
system must be substituted by 

9^ = /* . r {z., ~ x.,f exp(r,^^4-^) exp(-r,^^)0,dx. (17) 
6'i=exp(-/, J / exp(r. )exp(-ri— )dxi) (18) 

'Htot i J —oo i S* i 

Here i ^ 1. In order to simplify this system it is necessary to take into account 
the following statement 

• The size spectrum of the droplets of those sorts for which the formation 
of the droplets has been already described can be treated as the set of the 
monodisperse spectrums with some coordinates (depended on the sort of 
the centers). 

The solution of the system of the condensation equations was necessary to get 
the number of the droplets and the coordinate of the spectrum. 

Then the system of the condensation equations can be reformulated in the 
following way 



Ntot 1 „3 



= / (z, - a:,)^exp(r, ^ ' ' )exp(-r, .^s^^)e,dx, (19) 

9i = exp(-/, i / exp(ri )exp(-ri )dxi) 

Vtot i J — oo S+ i i 

(20) 

Here i ^ 1, Ntot i is the total number of the droplets formed on the heteroge- 
neous centers of the first sort and zi i is the coordinate of spectrum of the sizes 
of the droplets formed on the heterogeneous centers of the first sort. 
Actually 

Ni tot ~ tot 

If the last equality isn't valid it means that the supersaturation is going to attain 
soon the maximum and then to fall. So, no other formation of the droplets can 
be observed. 
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One must solve the last system for all sorts of the heterogeneous centers and 
construct the sequence {C* i}i>2- Ordinary the sequence will be the old one, 
but the concrete values are changed. 

Again we can use only the lowest C* i, i-S- the second sort. Then we must 
repeat this procedure again and again. The question how to solve these systems 
will be discussed later. 

After condensation on jq sorts of the heterogeneous centers has been already 
described and we are going to describe condensation on the centers of the sort 
i the following system must be solved 



= . / (z. - x.,f exp(r. '^(^^ ^* ' ) exp(-r. ^ ' ^ ' )e,dx. 



Ntat j ^3 

I 



= eM-f* / exp(r, "^"; exp(-r. "' ' ^'f' "°° '^'' )dx,) 

(22) 

Here i > ig, Nt^t j is the total number of the droplets formed on the centers of 
the sort j (it is near to r]j tot), Zj i is the coordinate of these droplets. 

System ( pl| ) - ( p^ ) is rather complex because the values Zj i are unknown 
ones. The next statement simplifies the problem and reduces it to the two more 
simple problems: 



Wtot 3 -,3 



(21) 



• System (21) - ( p2D can be with a high accuracy reduced to the following 

system 

5. = /* ^ r {z^ - x,)3exp(r,^^i^l— ^)exp(-r,-^)0,da;, (23) 
9. = eM-f* r exp(r,^^fc^) exp(-r,-|^)dxO (24) 

"^tot i J —oo i i 

where the function H satisfies the following equation 

^^n+Y,^^{z + 5z,f (25) 

(26) 

a. at 

and Szj are some fixed and already known distances between the monodis- 
perse spectrums. 

Here the following statement is also taken into account 
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• The following approximation during the period of the intensive formation 
of the droplets on the heterogeneous centers of the given sort (PIFDGS) 
is valid for the nucleation rates Ji 

Jz = Jz(r/tot,,f^*)exp(r,^^^^^)^ (27) 

where 

r. = -^^.^lc^o. (28) 

and the index "*" corresponds to i. 

The validity of the last system of equations is based on the simple fact that 
the rate of the vapor consumption is violated by the relatively essential deviation 
of the supersaturation and the rate of nucleation is violated by the relatively 
small violation of the supersaturation corresponding to the violation of ^ . 

We shall call 17 the external supersaturation. Our next task is to determine 

it. 



4 Iterations for the external supersaturation 



System ( |25| ) - (^6|) forms the ordinary differential equation (of the Abel's type) 
and doesn't allow the analytical solution. The initial condition is the following 

where tmit is the time of formation of the last known peak of the droplets 
distribution. 

Meanwhile one can construct the iterations on the base of the small param- 
eter 

^_ T.,<^oNJtot[z + 5zJf 

We suppose that 

dt^ - 

and can announce the following statements 
• The iterations defined according to 



(29) 



r ■r-^ Nj tot , , J- ^3 



converge to a unique solution (an ordinary property of a differential equa- 
tion)! 



^We shall put for simplicity z^qj = 0. 
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• We need to investigate the region where 

dn „ 

— > 

dt ~ 

When ^ = at the final of the period under the consideration we have 

ActuaUy we can remove Unit to the time of formation of the first peak and 
consider all 6zj as the positive values. Then 

d\maxn(^k) - mQxl7(fc+i)| ^ ^ 
dSzj ~ 

The previous inequality in such a reformulation can be also justified. An- 
other inequality can be also taken into account 

Certainly, tinu can not be put earlier than the moment when <f> = 0. 
So in the worst situation the form of the differential equation is the following: 

At= — - + —^z^ (30) 
a dt Uoo 

where A is some constant and Ntot is some value. 

If $ is the polinom then all iterations are calculated analytically and have 
the form of polinoms. Their convergence is rather high. The value of a relative 
error in the maximum of Q 

^ ^ \maxtfl(^i) — maxtfl\ 
' ~ maxtfl 

is less than 

\ i+l 



\maxt^(i) - maxtfl^i+i)] 



i+l) 



maxt^( 

and can be easily calculated on the base of (^0[). As a result one can get 
A2 < 4.2 * 10^2^ A3 < 8.68 * 10^3^ 

Note that the requirement for the choice of the number of iterations is the 
following 



11 



Really, Ti is the big parameter of the theory and goes to infinity, but actually 
Ti ~ 50 which is necessary for the intensive formation of the droplets (/* ~ 
exp(— consfF^^^)). So, we need to calculate only the few first iterations. 

Figure 1 shows the form of the iterations. Note that the evident renormal- 
ization can lead to the differential equation without parameters 

with the previous initial condition. 

The forms of the iterations /i - are drawn in Figure 1 where the inex points 
the number of the iteration. Note, that only the behavior until the maximum is 
interesting for us. It can be seen that can not already be separated from f^. 

One can also introduce another approximation which is more accurate. This 
is the approximation of the ideal surface. In this approximation the term z'^ is 
split between the surface term and the linear term z. Return to the previous 
scales of the variables. The iterations are constructed by the following way 



T C?Z(fc+i) Nj tot , . 7 \2f , f \ 

= -dt + ^ + -^^^O (^(fe+1) + ^^3) 

'i<io ° 



The rate of convergence of such iterations is higher than of the old ones, but 
there is no opportunity to calculate iterations analytically. Hence, this way of 
the approximation can be effectively used only at the last step. 

Now we are going to solve system ( |23|) - ( p^ ) . This system is more complex 
than the system already considered in pi because the behavior of Vl during 
PIFDGS is essentially nonlinear. One can not linearize Vl as the function of time 
and X during PIFDGS. But some important conclusions about the behavior of 17 
can be made. Consider an arbitrary sort of the heterogeneous centers i on which 
the droplets have been already formed. Introduce the characteristic width Aa; 
of the peak which is going to be formed. Certainly, 

Ax < z, (32) 

The value of gi near z ^ allows to give the estimate 

Ni tot 3 
9i 

We can also give the analogous estimates for the derivatives 



dgi ^ 


Ni tot 2 


dz 


Zr 


(Pgi 


tot 


dz^ 


Zi 

^00 
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d^9i ^ Nj tot 

The actions of this terms on 17 are given by the following terms (i.e. the corre- 
sponding terms in the Tailor's series) 

dgi Ni tot 2 A 
• z. Ax 



dz 

d^g, N, 



-z^iAxY' 



dz^ 

Due to (52|) one can see that 



' '"'(Ax 



d^gi Ni totf^_^3 



n„ 



f^c. 

The action of ^ can be cancelled by ^ because they have the different signs. 

So the action of must be taken into account. The action of can be 
neglected. 

As the result we justify the following approximation for the behavior of Q 
during PIFDGS: 

Y 

— Q = cx + Ix"^ + const (331 

with the two known constants. 

The characteristic length of formation of the spectrum by means of the vapor 
consumption by the droplets in this very peak is given by the following formula 

The halfwidth of the spectrum initiated by the external supersaturation is given 

by 



i I (^^2 I Q — Qjnax I 

when formation of the droplets occurs near the maximum of Q. 

When Ax ^ Xp we can use the standard iteration procedure described in |^ . 
When Ax S> Xp or Ax ~ Xp we must use the modified method of the steepens 
descent [|| which is described in the next section (it covers all situations). The 
case Ax 3> Xp is extremely simple (C coincides with fl). 
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5 Modified method of the steepens descent 



In this section we shall forget about the several types of the heterogeneous 
centers and about the several peaks in the spectrum of the droplets sizes. We 
shall consider the model situation when we have the nonlinear external source 
of the supersaturation of the square form and the process of the homogeneous 
(or the heterogeneous) condensation. We shall start with the homogeneous 
condensation and then the generalization on the heterogeneous case will be 
given. According to the results of the previous section this model covers all 
possible situations. 



5.1 Homogeneous condensation 

In order to obtain the form of the function (^{x) the balance equation must be 
introduced. For the homogeneous condensation one can obtain the following 
balance equation 

^ = C + 9 9{z) = ^^2^ f dx{z - xff{x) (34) 

The value of can be chosen as the moment of the maximum of the supersat- 
uration. 

The term *y* is approximately equal to 1 and can be missed. 

The value of the ideal supersaturation can be presented in the following form 

$ = + $*(c2; + Za;2)/r (35) 

with the two parameters 

as it was done in The negative value of / corresponds to the nonlinear 
character of the external conditions. The form of the ideal supersaturation ( ^5|) 
can be treated as the Tailor's series cut of on the second term. The result of 
[|j states the validity of (35) in the case when the nonlinearity of the external 
source is induced by the vapor consumption of the already formed spectrums. 

The two assumptions justified analytically can be introduced. The first one 
establishes that 

• the leading role in the vapor consumption belongs to the super-critical 
embryos, i.e. to the droplets. 

The second assumption states 

• the quasistationary state of the embryos in the near-critical region. 
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It must be valid during the period of the essential formation of the droplets. 
Let be the characteristic time of the variation of the stationary state in the 
near-critical region and let ts be the time of the relaxation to the stationary 
state in the near-critical region. Then the required assumption can be written 
as the following one 

<C > ts (37) 

These assumptions will be valid in all situations of the homogeneous case and 
in the heterogeneous case except the situation when almost all heterogeneous 
centers are exhausted at the final of the period of the essential formation. In 
this situation the result of the period of the essential formation of the droplets 
is obvious: the number of the droplets coincides with the total number of the 
heterogeneous centers. The form of the spectrum is unessential during the period 
of the essential formation of the droplets and the spectrum is the monodisperse 
one after the supersaturation begins to fall. 

Note that the second statement may be not valid in some situations when 
the supersaturation doesn't depend of the process of the vapor consumption by 
the droplets (it is negligible). As it can be shown in these situations the effects 
of the nonstationarity aren't important for the results of the process^. 

Now let us analyze the behavior of g{z). At first we shall analyze it only 
qualitatively. We shall extract the three regions. Let us put / '--^ 0. At first 
we shall define the period of the creation of the main consumers of the vapor 
during the period of the essential formation of the droplets. Then the behavior 
of g in the first iteration can be seen from 

/•oo 

g{z)r^exp{z) y^exp{-y)dy (38) 
Jo 

One can note the following important fact: 

• The subintegral function in ( p8| ) is sufficient only when 

1 ~ 2/mm <y< Vrnax 8 y = Cp (39) 

This region is the region of the creation of the main consumers of the vapor 
during the period of the essential formation of the droplets. It is moving 
in time along the p-axis (and the i-axis) with the velocity of the droplets 
growth^. When the last inequalities aren't valid the subintegral function 
is negligible. The period which corresponds to (|39| ) plays the main role in 
formation of the spectrum. 

Now we shall extract the other regions 

^This case requires the special consideration which can be done and shows that the effects 
of the nonstationarity are cancelled. 

■'Only until the end of the process of the creation of the droplets because the deviation of 
the supersaturation from the ideal one essentially transforms the definitions of the boundaries 
of this region 
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• During the period of the essential formation of the droplets one can extract 
the "initial" region according to 

- 9 < cz < czb (40) 

The value of Zf, must be established with the help of the additional condi- 
tion 

« 1 (41) 

where / is the size distribution of the droplets and /i is the imaginary 
distribution formed at the supersaturation $ instead of C- The values 
corresponding to this region will be marked by the subscript " in" . 

• The estimates for g{z) give the above estimate for the duration of the pe- 
riod of the essential formation of the droplets. The period of the intensive 
formation of the droplets finishes when cz ^ 1. It starts when cz ~ — 1. 
One can define the peak of the period of the intensive formation of the 
droplets according to 



— c 



< z < c-^ (42) 

The successful solution of the system of the condensation equations on the 
second step of iterations shows that the subintegral function lies essentially 
inside the initial period during the intensive formation of the droplets. 
Then one can conclude thalj^ 

- 2/mm < CZb (43) 

If / < (namely this case must be considered), then the hierarchy becomes 
even more strong. 

So, the region of the creation of the main consumers of the vapor during the 
period of the essential formation ofthe droplets belongs to the initial region. 

5.2 Approximation of the spectrum 

One need to construct some approximation for f{x) at the period of the essential 
formation of the droplets. The problem can be reduced to the suitable expression 
for C ~ C* in standard exponential approximation 

fix) = /, exp(^^i^^) exp( ^(^*^-^*) ) (44) 



^and even 

'~ymin -f- 1 ~ 
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One can decompose C ~ C* into the Tailor's series in the neighborhood of 
z = 0. The expressions for the derivatives are the following ones: 

72a ^2, 



-6 f J{x){z-x)dx (46) 

J — OO 



The value d'^g/dz^ is proportional to the number of the already formed 
droplets. This value is the main result of the process of condensation and the 
subject of interest of the theoretical description. According to the iteration 
procedure one can note that 

• 1). At the first step of the iteration procedure the expression for the total 
number of the droplets isn't correct at all. The relative error goes to oo. 

• 2). The second step of the iteration procedure is the final one and it gives 
the almost precise expressions for all main parameters of the process of 
condensation. 

Hence, one can introduce the procedure without the account of (|4^), (^8|). 
Then one can can restrict the Tailor's series by the first two terms. This restric- 
tion leads to the following approximation 

fix) = f*eMr^ U=o ^J^M^;^ — ) = /™^-p(r^l-o^) 

(49) 

This approximation is valid only inside the peak of the period of the essential 
formation of the droplets. It fails outside this region. The standard method 
of the steepens descent spreads approximation ( ^ ) on the whole period of the 
essential formation of the droplets (including the " initial" region and the region 
of the creation of the main consumers of the vapor during the period of the 
essential formation of the droplets) which leads to the essential relative error. 
Moreover, approximation ( p9[ ) isn't necessary because in the initial region the 
correct expression for the spectrum has been already established 

/(x) = /,exp(r ^('^^~^* ) (50) 
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5.3 Equations on spectrum parameters 



To obtain the system of equations on the parameters of the spectrum one can 
differentiate the balance equation for the condensated substance at 2 = 0. Then 
the foUowing equations appear: 

$(0)=5m(0)+gex(0)+C(0) (51) 



d$ dg,n , dgex dC, 



The subscripts " in" and " ex" mark the part of g referred to the droplets formed 
inside the "initial" region and outside it. Equations ( |50| ) and ( ^5|) lead to 

3 

g^n = z^-'ui (54) 



i=0 

2 



^=^z2-»(3-j)a, (55) 



i=0 



1 



= ^zi-X3-«)(2-i)a, (56) 



dz2 

In these equations the values ai are defined according to 

= -i.o ' -m /* / 2;* exp(ca; + Zx^)cia; (57) 

To obtain the algebraic system of equations one can substitute the integral term 
in the last expression by the one of the Boyd's estimates 

' < exp z^) / exp(-r)dt < ' = z>0 

V^^T^ + (tT - 1)Z - ^ 'J, ^ ^ - ^(tT - 2)z2 + TT + 2z 

(58) 

Analogous procedure can be proposed for the values ffex , ^^f|^ , ^jff^ • Ac- 
cording to ( p9| ) the expressions for gf,x{Q), dg^x/dz \z=o-, d^gex/dz^ \z=o have the 
following system 

fO 

gex{0) - fm / expi-{ — )^)x^dx = /^X^Ag (59) 



dge 



dz 



=o=3/,„/ eM~{ — f)x^dx = ?,fmxl\2 (60) 

J zt Xp 
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--f^ U=o= 6/m f exp{~{ — )^)xdx = e/mXpAi (61) 
az J Xp 

where 

reexp{~e)d^ eo = - (62) 

and Xp is some parameter like the characteristic halfwidth. 

As the resuh one can see the system of the three algebraic equations which 
can be solved by the standard numerical methods. The solution of this equation 
gives all main characteristics of the process. 



5.4 Linear source 



In the case I = the final expressions become more simple. For Zb one can take 
here the value ~Xp. Then one can obtain 

This induces the following expressions 



_3! 

J=0 



Sm(0) =a3 = -/*c-4^^^— yy(-l)J"(-cXp)3-J"exp(-ca;p) (64) 



U=o= "2 = /*3c ^^{-^y ^ e^p{-cxp) (65) 

^ U=o^ 2ai = -/*6c-2 Y^i-^Yi-cXp^^' eM-cXp) (66) 
The values of become the universal constants 

Aj = y exp(-x^)x'da; (67) 
Introduction of the connection between and Xp 



c?z^ Fxp 
leads to the transformation from (|6^) - (|66|) to the following form 



(68) 
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^ = /,c-33^ -^(~l)^(-cxp)2-^ exp(-cxp) + 3fraxlX2 (70) 

= -/*c-26^(-l)^"( 
From (|o|) it follows that 

/* = %^[c-33^7;j^(-l)-'(-ca;p)2-^"exp(-cXp) + 3vI/:r3A2]-i (72) 



_ $.c{c-^6Ej^o(-l)J(-ca;p)i~Jexp(-cXp) + 6xgAi^} _ 2$, 
" r{c-33E -=0 (2%(-l)^(-c^p)'-^' exp(-cxp) + 3VI/X3A2} " 

To solve the last equation one can simplify it by the following transformation: 
Consider as some known value. It can be found from (|69|). In (|6^), (^ in 
some rough approximation the last term is negligible in comparison with the sec- 
ond term and one can spread the approximation of the "ideal" supersaturation 
for all this region. Then we immediately get 

* = exp(-l) 

Another simplification is to notice that at the characteristic scale cxp ^ 1 the 
value of exp(— ca:p) is near to exp(— 1) and can be linearized. When and T 
are the given values then equation ( [73[ ) for Xp is the ordinary algebraic equation 
of the power 4 and can be analutically solved. Then ( |7^ ) gives the value for 
fm and, hence, for and F. This completes the current step of the iteration 
procedure for The structure of solution in the nonlinear case is the same 
one. 

System (|69|)-(|7ll) allows some modifications which are available for the non- 
linear case also. Note that the subintegral function f{x)x in the expression for 



2$ 

-cXpY^^ exp(-cxp) + 6/„XpAi - 



(71) 



where 



and (f7l[) transforms to 
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(Pg/dP is essential only in the finite region of the period of the essential for- 
mation of the droplets. The known approximation ( p9| ) also leads to the same 
conclusion. Then approximation ( p9| ) can be formally spread on the whole re- 
gion of the essential values of the function f{x)x. Then equation ( |7l| ) leads to 
the following expression for Xp 

In the nonlinear case equation (^l|) has the following form 

^ = ixlU - 2^ (75) 

This bisquare equation gives Xp dependence on f^. Then the second equation 
of system - (|7^) becomes the close equation on the supersaturation and 
can be solved by iterations. These iterations can utilize the sharp dependence 
of exp(— i^c) on the supersaturation. The zero approximation can be chosen as 
T ~ 1, $* = 1. Then the second approximation gives rather a precise result for 
all parameters of the spectrum. Certainly, the amplitude fm must be obtained 
by ( |7l| ) instead of the explicit expression from the classical theory of nucleation 
on the base of the supersaturation. 

The system of the condensation equations can be even more simplified if one 
notices that approximation (|50| ) can be spread over the whole region of integra- 
tion in the first and in the second equations of the system. This simplification is 
valid as far as the subintegral functions in the first and in the second equations 
are essential only in the " initial" region of the period of the essential formation 
of droplets. Finally, the system of the condensation equations has the following 
form 

—;—\z={)—'if* f x^ exp(cx + lx^)dx (76) 

= 6Aix2/„ - 2^ Ai = y" yexp(-y^)dy = 0.5 (77) 
One must take into account the difference between fm and /* which is given 

by 

= cxp( — \z=o / x^ eyip(cx+lx'^)dx/ / x'^ exp(cx+lx'^)dx) (78) 

It is essential that equations ( |76| ) and ( |77| ) are the separate ones. Equation^ 
^ can be solved by iterationfj^ The total number of droplets can be obtained 



as 



^The integral can be calculated explicitly and we have the ordinary algebraic equation. 
^As the algebraic equation due to the sharp dependence of /, on the supersaturation 
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N = frnVT^XpHoo (79) 

The value of parameter I is negative. But one can prove that the relative 
error of the last expression 6 increases when I decreases in the absolute value 

^1^1 >0 (80) 



dl 



In the linear case 



5(0) = ^ (81) 

Xp = c-i(4exp(l))i/4 (82) 

N ^ -^^^^(4exp(l))V-V^ 

6rexp(l) °° ^ ' 

Then N is 1.19 times greater than the result of iteration procedure at the 
second step and 1.03 times greater than the result at the third step. It can 
be regarded as the practically precise result. All mentioned above estimates 
obtained for / — remain valid in the nonlinear case. 

Figure 2 illustrates these constructions in the linear case. This picture shows 
the contributions of the various forms of the spectrums into the values of 5, 
dg/dz, d^g/dz^, d^g/dz^ (from upper till lower axes of the coordinates) atQ 
t = t^. Here the subingral expressions obtained in the various approaches are 
drawn. The upper picture shows the subintegral functions in the expression for 
g a.t z — 0, i.e. f{x). The second picture shows the subintegral function in the 
expression for dg/3dz at z = 0, i.e. f{x). The third picture shows the behavior 
of the subintegral expression for d^g/6dz^ at 2; = 0, i.e. xf{x). The lower picture 
shows the subintegral expression for the number of the droplets or for d^g/6dz^, 
i.e. f{x). For f{x) several different approaches are accepted. The curve "a" 
corresponds to the precise universal solution, the curve "b" corresponds to the 
spectrum calculated on the base of the ideal supersaturation (i.e. the first 
iteration), the curve "c" corresponds to the gaussian form of the spectrum and 
the curve " d" can be drawn on the base of the spectrum of the gaussian type 
with the cube correction term (see later). 

It can be seen that for g the suitable approximations can be given by the 
first iteration. For dg/dz the suitable approximations can be given by the first 
iteration and the gaussian spectrum with a correction term (with the appropri- 
ate cut-off). For d^g/dz^ and for d^g/dz^ the suitable results can be obtained 
from all approximations except the first iteration. 



*In the homogeneous condensation the choice of i* corresponds to the maximum of the 
supersaturation. 
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The adopted approximation has some resources of the modification inside 
itself. Namely, it can give the value of the correction term and the value of the 
boundary of the region of the applicability of this approximation. Really, let us 
add the next term to the square approximation. Then we have 



Xp Xl 



with some parameter xi. The value of xi will be established from the fact that 
the model curve for the supersaturation 



Cappr — C* p( (') (^)) 
i Xp Xl 



must only touch the line of the ideal supersaturation at some point, which will 
mark the boundary between the initial region and the region of the applicabily 
of above given approximation. 

Figure 3 illustrates this construction. Figure 3 demonstrates the form of 
the modified spectrum with the cubic term taken into account. The upper 
picture shows the behavior of the normalized deviation of the supersaturation 
as function of z for the different models. The value Sid corresponds to the ideal 
supersaturation, the value Ssq corresponds to the gaussian approximation, the 
value Smod comes from the consideration of the gaussian with a correction term. 
The curve Smod touches the straight line Sid- In the lower picture the gaussian 
spectrum f gauss and the modified spectrum fmod are drawn. The scale of x is 
chosen to have c — I. It is seen that the difference between these two spectrums 
is rather small. 

In the case of the linear source it follows from the condition of equality of 
the derivatives of Cappr and $ that 

-I - {-\-~pW - 

After the substitution of into the condition^ Cappr{xb) = ^{xb) we shall come 
to 

-2 + (4-12Xp2)i/2 



Xb 



for the boundary between the initial and the central regions. 

By the same way one can get the number of the molecules in the droplets by 
the separate integration and, thus, fulfil the next step of the iteration approxi- 
mation. 



'We assume that approximately 't = exp(— 1). 
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5.5 Heterogeneous condensation 

In the case of the heterogeneous condensation one can obtain the foUowing 
system at the period of the essential formation of the droplets 

^ = ( + 9 9{z)= f dx{z-xff.9exp C^^^*'^^ ) (84) 
fc* = exp( / dxf^exm ^ )) (85) 

mot J-oo 

The straight development of the method utilized in the homogeneous case 
leads to the construction of the square approximation for the function 



If t r ^ f , -r(<i>,-c(x)) ,, -r{<f,-c{z)) , 

V' = /*exp / dxf^expi ) exp ) 86) 

which comes from the substitution of the second equation into the first one. 
The square approximation for isn't suitable because the factor 

exp(-— / dxf^ exp( ^^^* — ^)) 

Vtot J-oo 

in the last expression has a very complex behavior. This term goes from the 
expression for 9. This allows to introduce the following method. Notice that 
the iteration procedure in the investigation of the heterogeneous condensation 
gives the appropriate result on the second step The second approximation 
is obtained on the base of approximation which doesn't contain any functional 
dependence of 6 upon z. Then one can solve the pscudo homogeneous equation, 
i.e. equation ( p4| ) with 9=1 and then calculate 9 according to the final step 
of the iteration procedure |^ . The pseudo homogeneous equation will be solved 
by the method of the previous subsections. This solution will give the main 
parameters fm^Xp. The function 6 is approximately given by the following 
expression 

9 = exp(-e(z6 - z)(t)' - e(z - Zb)0") (87) 



where 



/ exp(cx + /x^)da; (88) 



0" = /,!^ / exp(c2; + lx'')dx + U— I ~ exp{-{ — f)dx (89) 

Vtot J-oo Vtot J zb Xp 

This expression can be taken in the elementary functions due to (|58|). 

The number of the droplets which appears from in a very simple form 

N{z) = r]tot{l-e{z)) (90) 
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based on the conservation law for the heterogeneous centers. 

One can introduce another approximation. In the final formulas one can 
spread approximation (^9|) over the whole period of the essential formation of 
the droplets. Then the value of final can be presented in the following form 

dfinal = exp[-^^fmVT^Xp] (91) 
Vtot 

and the total number of droplets is given by 

71 

Ntot = Vtot{^- exp [ — frn V^Xp] ) (92) 

Vtot 

The error of the above presented method is rather small, the accuracy is 
practically the same as the accuracy of the iteration procedure |^ . The advan- 
tage of this method lies in the fact that the iterations can't be calculated in the 
case of the nonlinear external conditions. 

Note that in the above presented procedure the cross influence of the ex- 
haustion of the heterogeneous centers can not be observed. Really, suppose 
that there are many heterogeneous centers of the first sort and a few centers 
of the other sort. Certainly, the second sort can not act upon condensation 
on the first sort in reality. But in final formulas the second sort acts as there 
is no exhaustion of the second sort centers. The values <&' and $" are based 
on some model for supersaturation. Only the parameter Xp is determined by 
the consumption of the vapor by the droplets. But it is calculated on the base 
of the total values of the heterogeneous centers. So, we must reconsider this 
procedure. 

The process of exhaustion of the heterogeneous centers makes the inequalities 
for the hierarchy of the regions even more strong. So, we can repeat the above 
presented procedure with the account of the exhaustion of the heterogeneous 
centers. 

For the approximation of dj {j marks the sort of the centers) in the " initial" 
period one has 

r f* j 

tij = exp explcx) 

c rjtot j 

Then this expression must be introduced into all expressions for a^. These 
expressions now are the following ones 

ai J = (-1)' ' J ^j I exp(ca; -I- ^x^) exp[ — — 2^ exp (ex)] 

i!(3-i)! c rjtot j 

When the last expression is going to be calculated we shall also use the 
steepens descent method and decompose 9 into the Tailor's series near x = Zh. 
In the zero approximation 

ttj 7 = (-1)*-T77^-— tt/*7 / exp(ca: + lx^)dx exp[- exp(czfc)] 
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has the same analytical structure as the previous one. Also one can decompose 
exp(ca;) near exp(czb) and restrict this decomposition by the zero and the first 
two terms of the Tailor's serials. 

The further generalization is rather simple. We must include all times 
when /» j appeared. 

In the case of the continuous spectrum we must substitute by J dw where 
w is the activity of the heterogeneous center. 

Note that the effect of the cross influence of the exhaustion of the hetero- 
geneous centers doesn't play any essential role in the dynamic conditions. This 
effect isn't specific for dynamic conditions but takes place in the situation of 
the decay type. The reason is that in the situation of the decay there are 
many characteristic widths at one and the same (initial) moment of time for the 
given initial value of supersaturation. In dynamic conditions all halfwidths are 
approximately equal to c~^. 

The alternative opportunity is to consider exp(/a;^) as the correction term 
and to keep exp(— exp(ca;)) under the sign of integral. 

One can say that in the situation of the dynamic conditions the process of 
the exhaustion (the probability to be exhausted) occurs only due to the activity 
of this very sort of the nucleus (at the given supersaturation). 

6 Concluding remarks 

The applicability of the quasistationary distribution as the boundary condition 
for the kinetic equation is based on the fact that the characteristic time of the 
establishing of the stationary distribution in the region i/^^'^ < 3^1^^ is negligible 
in comparison with the characteristic time of the variation of the supersatura- 
tion up to the value C/^i- This statement is valid in all situations except the 
situations when there occurs the total exhaustion of the free heterogeneous cen- 
ters and the result is obvious - all centers of the given sort became the centers of 
the droplets. This fact can be proved analytically. Also one can prove that the 
main consumers of the vapor are the super-critical embryos, i.e. the droplets. 

The quasistationarity may not be valid when Ax ^ Xp. One can show that 
in this case the effects of nonstationarity are cancelled. 

Consider how to include situation (^3|) in the general consideration. Note 
that in the iteration procedure |6j as far as in the modified method of the 
steepens descent the structure of the calculations has one remarkable feature: 
the final expressions for 9 and N are given on the base of the approximation 
obtained with the help of the "ideal" laws. In these "ideal" laws the analytical 
structure of the expressions for the different sorts is one and the same. It turns 
to manifest some features specific to the sort of the heterogeneous centers (due to 
their number) only at the final iterations. So, the generalization is rather formal 
and simple. It is necessary to have one common system of the coordinates a;, 
z and to calculate directly the iterations (or to use the modified method of the 
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steepens descent). For 0(2) i we get the separate expressions: 



Here no cross influence of the exhaustion of the heterogeneous centers of the 
different sorts can be observed. One must use the modified method of the 
steepens descent with the cross influence of the exhaustion of the heterogeneous 
centers as had been made for the case of the decay. The result of condensation 
on the spectrum of activities will show that the form of the spectrum is near 
the universal one. So, one can use it as some ansatz and formulate the algebraic 
equations on the parameters. 

The account of the nonisothermal effects is rather simple. The intensity of 
the droplets formation is changed in ln(C + 1) times (i.e. in the approxi- 
mately constant number of times during the period of the intensive formation 
of the droplets). So, the description of the process can be attained if we assume 
that the consumption of one molecule effectively leads to the consumption of 
(^2 + 1) molecules, where fc2 is some known parameter So the description 
is reduced to the obvious renormalizations. 

The global description of the evolution is more complex. One can show that 
the already known procedure of the iteration approximations for the external 
supersaturation (where z in all terms except ^ is supposed to be known) is 
successful^ in this case also. The application of this iteration procedure allows 
to cancel the restriction ^2 ^ 1 used in 

Practically one needn't to fulfill such a complex procedure as to solve many 
system of the integral equations many times. Ordinary the following procedure 
must by accepted: 

• Put C, — ^ and determine the moment ti of the exhaustion of the different 
sorts of the heterogeneous centers. 

• Establish the maximum of under the monodisperse approximation with 
the already defined coordinates of the peaks and iVj = rjtot i 

• Reconsider formation of the droplets of that sorts of the heterogeneous 
centers for which 

z{ti) > z(t(max fl)) — Ax 
by the modified steepens descent methodpl 



^"Here the initial approximation corresponds to the ideal law of evolution 
^^Note that formation of the droplets on the sorts for which 

z{ti) > zitimax n)) + Ax 

can be treated in the pseudo homogeneous manner 
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As far as the moments of formation of the droplets of these sorts are near 
t{max ft) the action of formation can be presented by the direct summation of 
all /* i for these sorts in the expression for g and its derivatives. The last step 
solves the problem. 
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Figure 2 

Subintegral functions in different approximations. 
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Modified spectrum with cubic term. 
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